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بکار گیری روش المان‌محدود در تحلیل همسته سدهای خاکی و برآورد فشار آب‌منفذی * 


فرهود کلدی,(۱) فریده حسین نواد!۳؟ 


چکیده افزایش قابل ملاحظه فشارها یآب منفذی هسته حین ساعت سدهای نعاکی می‌تواند منجر به بروز پدیده‌ی شکست هیدرولیکی 
هنگام آبگیری سد گردد. لذا اطلاع کامل از الگوی ایجاد اضافه فشا رآب حفره‌ای در هسته سد امری ضروری می‌باشد. در تحفیق حاضر با 
بکارگیری روش المان محدود در تحلیل همبسته سدهای نعاکی» مقادیر فشا رآب مفذی و تغییرمکان‌های ایجاد شده حین سانخحت سد درونگر 
با استفاده از کد نکاشتته شده در محیط فرترن برآورد گردیده و نتایج با داده های ابزاردقیق مقایسه شده‌اند. درنهایت با اسفاده ا زآنالیز 


با ز؟ گشتی مقدار واقعی نفوذپذیری هسته ونسبت فوذپذیری افقی به قاثم مشخحص شده است. 


واژه های کلیدی سد خاکی, معادلات بیوت تحلیل کوپل فشار آب منفذی. 
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* تاریخ دریافت مقاله ۹۵/۱/۲۰ و تاریخ پذیرش آن ۹۱/۲/۲ می‌باشد. 
(۱) استادیار دانشکده مهندسی عمران- دانشگاه تبریز. 


(۲) نويسندهٌ مسئول. دانشجوی دکتری. مهندسی عمران- سازه های هیدرولیکی. دانشگاه تبریز. ‏ عطمه.0مطهو۵06۵زممصزهومط_طعننط باتفصظ 


۳ بکارگیری روش المان‌محدود در تحلیل همبسته سدهای نعاکی و برآورد فشا رآب‌مفذی 


مقدمه 
امنیت سازه های بزرگ. حفاظت زیست محیطی و 
کاهش خسارات وارده در بلایای طبیعی مستلزم فهم 
درستی از دلایل و مکانیزم‌های فرایند تغییر شکل سازه- 
هامت: [1 شهار اد یله مبانهای ممخاسی تخت 
دست بشر می‌باشند که به طور مستقیم با آب و نیروهای 
مخرب آن در تماس هستند. در یک دسته‌بندی کلی 
سدها به دونوع بتنی و خاکی تقسیم‌بندی می‌شوند. 
سدهای خاکی و سنگریز به واسطه مزیت‌های نسبی که 
در مقایسه با سایر انواع سدها دارند» دارای حایگاه 
ویژه‌ای بوده و بخش عمده‌ای از سدهای احداث شده 
در سطح جهان را (در حدود ۸۳ درصد) به خود 
احتصاص داده‌اند [۲]. اطلاعات آماری اشاره بر این 
دارد که سدهای خاکی بیش از سدهای بتنی در معرض 
تخریب هستند و بیش‌تر خرابی‌ها یا آسیب دیدگی‌های 
سدهای خاکی در حین ساخت و يا اولین آبگیری بوده 
است [3]. یکی از مشکلات اساسی در سدهای خاکی 
افزايش اضافه فشار آب حفره‌ای در هسته سد طی 
مراحل ساخت و آبگیری می‌باشد که این مسئله می‌تواند 
منجر به پیدایش ترک هیدرولیکی در سد شود. لذا 
توانایی پیش‌بینی گسترش و توزیع (زائل شدن) اضافه 
فشار آب حفره‌ای در بررسی رفتار چنین سازه‌هایی از 
اهمیت بالایی برخوردار است [41]. روش المان محدود 
ابزار قدرتمندی برای تحلیل و حل مسائل در ساخت 
سدهای خاکی می‌باشد چون می تواند تغییرمکانهای 
داخلی هسته و پوسته را محاسبه کرده و در نتیجه توزیع 
تنش و انتقال بار داحل مقاطع سد را بدست دهد [۵], 
محققان بسیاری از این روش برای بررسی تغییرشکل‌ها 
و تنش‌های سدهای خاکی در مرحله ساخت و اولین 
آبگیری استفاده کرده‌اند [5-12] و [۶] و[1] . اما در اکثر 
موارد در تحلیل‌ها از روش تنش کل بهره گرفته شده 
[6,8,9,10,12,13] و در صورت استفاده از روش تنش 
موثر نیز اضافه فشارهای حفره‌ای ایجاد شده به واسطه 
ضریب پیش فرضی از تنش فرض شده است که تا 
ای تتخافظه کا رنه اشکدر ان آترن رشیگنه: کی 
همبسته (کوپل) که مبنای آن روش تنش موثر می‌باشد 
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ری مامتان اش ات ار هن و 
بررسی فشار منفذی ایجاد شده در هسته سدهای خاکی 
نیز به دلیل اهمیت فراوانش مورد توجه محققین زیادی 
بوده است. پاکزاد و حسینی [۱۵] با استفاده از نتایج 
پیزومترهای به کار رفته در پی و بدنه سد کرخه 
وضعیت سد را پس از آبگیری بررسی کرده و نشان 
دادند با توجه به ضریب فشار منفذی ماکزیمم حین 
ساحت و آبگیری» بروز پدیده شکست هیدرولیکی در 
این سد غیرمحتمل است. همچنین حسینی و همکاران 
[۱0]» افزایش فشار منفذی در هستهٌ سد کرخه را هنگام 
ساحت و آبگیری همزمان. مورد ارزیابی قرار داده و 
نشان دادند نتایج تحلیل‌های عددی در دوران ساخت 
نسبت به محور سد تقارن دارند ولی حین آبگیری مقدار 
فشار منفذی بالادست از حالت تقارن خارج شده و 
افزایش می‌یابد. نیرومند و همکاران [۱۷]» نیز با استفاده 
از نتایج ابزار دقیق و آنالیز برگشتی انجام شده با نرم 
افزار 6۸2 عملکره نید کرخه را پروسی کرفند و نشان 
دادند به دلیل طولانی‌تر شدن مدت اجرا مقادیر فشار 
منفذی حدودا نصف مقادیر در نظر گرفته شده در 
طراحی هستند. حبشی و همکاران [۱۸] فشار آب 
منفذی ایجاد شده در هستهٌ سد خاکی دوستی را با 
استفاده از آمار اندازهگیری شده و نتایج تحلیلی بدست 
امده مورد بررسی قرار داده و با استفاده از تحلیل 
بازگشتی نسبت نفوذپذیری افقی به قائم مصالح هسته را 
برابر با ۲ بدست آوردند. ملکی و علوی‌فر [19] 
رفتارنگاری مسجد سلیمان را به همراه آنالیز عددی 
توسط کد نرم‌افزاری 4.0 ۸6 مورد بررسی قرار داده 
و پارامترهای ضریب فشار آب منفذی و نسبت قوس - 
شدگی را محاسبه کردند. سروش و اخترپورا ۲۰] مراحل 
ساخت و اولین آبگیری سد خاکی مسجد سلیمان را با 
نرم افزار ۳۱۵*1572 مدل کرده و فشارهای منفذی 
بدست آمده را با نتایجم حاصل از ابزار دقیق مقایسه 
کردند و با استفاده از آنالیز بازگشتی نسبت نفوذپذیری 
افقی به قائم مصالح یتست . او رونت با 
همچنین نشان دادند در زمان‌های اولیه خاکریزی بدلیل 
کم بودن سربار و اشباع نبودن کامل خاک مقادیر اندازه 


سال سی و یکم شمارة دو ۱۳۹۷ 


فرهود کلاته- فریده حسین نژاد 


گیری شده فشار آب منفذی کمتر از مقادیر محاسبه 
تور انس تیا اف ای گرا ها کی الا ای دی 
مقدار کمتر می‌شود زمردیان و پوچی [۱1] طی یک 
بررسی با استفاده از نرم‌افزار 0 رفتار سد 
خحاکی مسجد سلیمان در حین ساخحت و اولین آبگیری 
را به روش عددی و به صورت همبسته تحلیل کرده و 
نتایج حاصل برای فشار آب منفذی» ضریب فشار آب 
منفذی و تنش کل را با داده‌های ابزار دقیق مقایسه کرده 
و نشان دادند مقادیر محاسباتی و اندازه‌گیری شده برای 
فشار منفذی در ترازهای پایین‌تر مطابقت بیشتری دارند 
ولی در ترازهای بالاتر مقادیر اندازه‌گیری شده بیشتر از 
تقادیر امخاشانی ات سین با شیر یلها اکن با 
دلیل ازدیاد سربارهای ناشی از افزايش ضخامت 
عاکرفای گنها رت رین سا تفای مد 
افتد که نرخ افزایش آنها وابسته به نفوذپذیری مصالح 
می‌باشد. از طرفی سرعت بالای خاکریزی یکی از 
عوامل افزايش فشار آب منفذی داخل هسته به شمار 
می‌رود. در مکانیک خاک تغییرشکل حجمی خاک ناشی 
از تغییر تنش‌موثر با عبارت تحکیم بیان می‌شود و 
سرعت تحکیم در خاک اشباع با سرعت خروج آب از 
خاک تعبین می‌شود؛ بنابراین فرایند تغییرشکل و 
جریان‌حفره‌ای در محیط‌های خاکی با هم در تعامل‌اند. 
مسئله تحکیم در مسائل ژئوتکنیکی همچون پی‌هاء 
توئل‌ها و سدهای خاکی بسیار حاثر اهمیت است و از 
اینرو تحلیل المان محدود همبسته در مسائل ژئوتکنیکی 
بسیار مورد توجه بوده است. پدیده تحکیم اولین بار 
توسط ترزاقی درسال ۱۹۲۶ برای تحلیل یک‌بعدی 
محیط‌های متخلخل مطرح گردید و به طور وسیع در 
فتخاسیه تتصته تا کب کاه رف[ اف نها ترگ 
[22] با توسعه تئوری ترزاقی به سه بعد. برای اولین بار 
فرمول‌بندی تحکیم کوپل(همبسته) و حل تحلیلی آن را 
بر اساس رابطه خطی تنش-کرنش اسکلت خاک و 
رابطه خطی دارسی برای جریان حفره‌ای در مسائل شبه 
استاتیک(2110]-00261) ارائه داد و مدتی بعد معادلات را 
به حالت دینامیکی نیز توسعه داد. در ادامه کار او 
زینگویچ ساده‌سازی‌هایی در معادلات انجام داده و 
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۳۵ 


بیشتر به جزئیات معادلات پرداخته است [23-25]. در 
چند دهه اخیر نیز تحقیقات فراوانی در زمینه حل عددی 
عاالات یه ووامسظاهای شملگ میورت گرفه 
الا باه طیل سای نوی تخل فتانکی 
هام اک یه کار ره انا [ هن ی 
هیر بر مساو ای تفر و بیان کرن رو -دتتایگر 
معادلات همبسته بیوت ساخت مرحله‌ای بدنه سد 
خاکی درونگر در ده لایه مدل می‌شود و مقادیر اضافه 
فشار منفذی و نشست‌های محاسبه شده با داده‌های 
ابزاردقیق مقایسه می‌شوند و از نتایج حاصل مقدار 
مناسب نسبت ضریب نفوذپذیری افقی به نفوذپذیری 
قائی جهت انطباق بهتر نتایج با داده‌های ابزاردقیق تعیین 


می شود. 


معادلات حاکم 
معادلات حاکم بر رفتار محیط متخلخل اشباع که دارای 
راک هراق هت نانک کر 
ریب وتان سای اوه تتاون سرا تفر 
ی 


جرم سیال و مفهوم تنش موثر به صورت زير بیان می- 


شود [23,24] 
)0 0 < رم - [ررترن + رمع + نتم ره 
)0 9 9 ِ نابات بتامم د ۳ ۳۹ 8 
۲ 12 1( -1) 0 0 
فا ۳ 
(۳ 0 -< ۵ + 


تیمها وان انا سین تاشته و نان 
دهنده جهت ها( ولو ولو ) در سیستم مختصات 
کارتزین می‌باشند. وت تانسور تنش کل در نقطه 
( ولو ولو ر) و علامت "ز," نشان دهنده مشتق ضمنی 
نسبت به جهت ز می باند(چ). در واقع در حالت 
کل قعاله (۱) سل هسدنه مسسعل اس اون 
نیروها در جهات ولو وتو ی که در حالت کرنش 
مسطح تبدیل به دو معادله می‌شود. در معادله تعادل 
کلی. چگالی کل ترکیب(م) به صورت زیر محاسبه می 
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۳۹ بکارگیری روش المان‌محدود در تحلیل همبسته سدهای تحاکی و برآورد فشا رآب‌مفذی 


سود 
() وم( -1) + ,۲۵ 2۶ وم( -1) + ,70 < م 
که 7 نسبت تخلخل. ,رم چگالی آب و وم چگالی 
دانه‌های خاک می باشد. در روابط (۱)و(۲). :9 نیروی 
حجمی وارد بر واحد جرم است. در معادله دوم ۷۲ 
نشان دهنده نیروی درگ ویسکوز می باشد که از قانون 
تراوش دارسی به صورت ۷-61۳ بدست می‌آید. در 
این رابطه 6 نفوذپذیری دینامیکی مصالح غیر 
ایزوتروییک و :۷ سرعت سیال است که در حالت 
ایزوتروپیک بودن مصالح نفوذپذیری محیط متخلخل 
در کلیه جهات یکسان بوده و نفوذپذیری دینامیکی برابر 
می شود. از طرفی ۷۳ معادل زع گرادیان هد کل در 
جهت ر می باشد که 0 هد آب حفره‌ای کل برابر 
انرژی پتانسیل سیال حفره‌ای است و از رابطه زیر 
پلشتتاهین. یلم 
(ظ( نب ۳ (نجرون< - و < ظ 
2 جهتی است که جاذبه اثر می‌کند و تور 
نشان دهنده سطحی است که هد کل نسبت به آن اندازه 
گیری می شود (با صرفنظر از اثر گرانش). عباراتی که 
در معادلات تعادل مومنتوم کلی ترکیب جامد-سیال 
حفره‌ای (معادله(۱)) و تعادل مومنتوم سیال (معادله 
زان تاو تا انس سای سس : 
سیال نسبت به ذرات جامد و عبارت جابه‌جایی این 
شتاب می‌باشند. در معادله تعادل جرم سیال پا پیوستگی 
(معادله (۳)» دو عبارت آخر به ترتیب مربوط به تغیسر 
چگالی و تغییرات دمایی جرم سیال می‌باشند که در 
حالت کلی می‌توان از آن‌ها صرفنظر کرد [25]. در این 
مدشن عون الک درا سانن و1 سیون باتک 
دانه‌های خاک می‌باشد. با در نظر گرفتن ۶ به عنوان 
تراکم پذری سیال که با جر »6 و و تراکم‌پذری ذرات 
جامد که با < و) تعریف می شوند. تراکم‌پذیری 
ترکیبی فازهای سیال و جامد می تواند به صورت زیر 
تعریف شود. 
)1( ۵( - 1) +۷2 2 ون( -0) +4 02 < 6 
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که 0 ضریب تصحیح فشار سیال حفره‌ای موثر 
بر ذرات جامد می‌باشد و برای مصالح ایزوتروپیک به 
صورت "1-39 < » تعریف می شود. برای اغلب 
اک ها 1 < ۵ در نظر گرفته می شود. در این رابطه 
بوک نشان دهنده مدول بالک متوسط نمونه خاک می- 
باشد. با جایگذاری )و 4 در معادله (۳) و حذف 
مات یل مه برفی ماه فرق به عمو رف ری 
ساده می‌شود: 
)۳ 0 - 2 + رغه + روا 

در دی وان زو تفای بو رل 
مقظلییل با درنظر گرفتق انلوکش انسکلتاستدبا 
سیال حفره‌ای در دو شرایط استاتیکی و دینامیکی می- 
تاشتت هر انزم فعاول عتفشار اسف دای :ربب غیت 
نسبی جریان سیال نسبت به فاز جامد :۷ و جابجایی 
یسکات سا سل ادها فحیت ره امس با فان ون 
صورتی که در مسائل دینامیکی. فرکانس مولفه‌های 
شتاب تحریک پایه پایین باشد. همانند ارتعاشات زلزله, 
در معادلات (۱) و (۲) عباراتی که زیرشان خحط کشیده 
شده است(عبارات مربوط به شتاب سیال) قابل صرفنظر 
کرفن هی باکتتد |25 با متفه ایسن عباوات متفر 13 
می‌تواند از معادلات حذف گردد. بنابراین معادله حاکم 
ساده‌سازی شده که شامل دو متغیر زلا و 0 می‌باشد به 
صورت زیر در می‌آید: 
۸( 0 < ,رم - نام ۲ ز و6 

در مصالح الاستیک رابطه تتش-کرنش با قانون 
رفتاری هوک.ع9 < 6 بیان می‌گردد. از طرفی با فرض 
تغییرمکان‌های کوچک. کرنش با استفاده از رابطه 
2 با تغییرمکان‌های گرهی مرتبط می‌گردد. 
بطوریکه 1 ماتریس کرنش-تغییرمکان می‌باشد. با 
جایگذاری اين روابط تنها مجهولات معادله. تغییرمکان- 
ها و فشارحفره‌ای می‌باشند. همچنین با ترکیب معادلات 
(۲) و( و حذف جملات قابسل چشض‌پوشنی؛ دومین 
معادله حاکم به صورت زیرخواهد بود: 
(٩)‏ 0< + زرغه + ز[(رطم - رنابم + زظ)زنک1] 


معادلات (۸) و(٩)‏ در کنار هم فرمول بندی ۷-0 


سال سی و یکم. شماره دو ۱۳۹۷ 


فرهود کلاته- فریده حسین نژاد 


یا به عبارتی معادلات همبسته بیوت را تشکیل می دهند 
که باید به صورت همزمان حل شوند. برای حل این 
معادلات نیاز به تعیین شرایط اولیه (10) و شرایط مرزی 
(30) است. برای تعادل مومنتوم کل محیط بر روی 
خی از مرن 1 تبسن کل #متحصن مین .برد دن ای 
که برای ,آ. جابجایی !. داده می شود. برای فاز سیال 


حل عددی معادلات 

برای حل عددی معادلات از روش المان محدود 
استفاده می‌کنیم. اولین گام در حل عددی. گسسته‌سازی 
مکانی و زمانی معادلات است. برای گسسته‌سازی 
مکانی معادلات از روش باقیمانده‌های وزن‌دان روش 
گالرکین استفاده می‌شود. گسسته‌سازی مکانی شامل 
جایگزینی متغیرهای معادلات حاکم یعنی ‏ و 9 با توابع 
تقفریب مناسب می باشد. 

8 2 2۱۱ < ۲ و 2۱۲ ز۱۳۵ 2 > با 
(۰( 

که آ ر] مقادیر گرهی این متغیرهامی باشند. 
"و نیز به ترتیب توابع شکل مناسب برای 
جابجایی و فشارحفره‌ای هستند. به دلیل وجود مشتق 
مرتبه دوم جابجایی و مشتق مرتبه اول فشار در معادلات 
حاکم به ترتیب نیاز به پیوستگی مرتبه اول و مرتبه 
صفر توابع شکل مربوطه می‌باشد. در این بررسی از 
توابع شکل مرتبه دوم برای جابجایی و مرتبه اول برای 
قکار انناژه منود 


تاو نار و9 چتا ب۷ بلا] > 17 


اسر | ۲ 


 ]* 0 8‏ بر 
۳ 0... ۱۳ 0 
)۱۱( ,۱۶,۰ ,۱] - ۸ 


با جایگذاری روابط (۷) و (۸) در معادلات تعادل 


و پیوستگی, شکل گسسته شده معادلات به صورت زیر 


خواهد بود: 
)۱۲( ([) < (0[)5] - (تا)[م] + ()[۳] 
۱۳( (۶2) - (5[)8] + (5)[ + (ت) 0[7] 


سال سی و يکي شماره دو ۱۳۹۷ 


۳۷ 


که[۷1], ماتریس جرم برابر 00 ۳(م(0۳ ] < [۲1 
مک ماتریس سختی 0[]8[060]  ]]8["‏ [م] 
0 ماتریس کوپل (اندرکنش) 0۳[060] *[8] | < [0] 
ماتریس انتقالی 00[م8] 0 [:5]] جِ 

5 ماتریس‌تراکم‌پذیری ۱۳40( + )۱۷ ] - 


و ]و 62] بردارهای نیرو می‌یاشند که به صورت زیر 


۱ 
1 
1 
1 


محاسبه می‌شوند. 

۲ )0 | +هه ۸( | ()۶ 
۷ 

۲ ۳(۲۵) + (06(ط,م,5) ۷۳ ۱۱۳(۲) 1 - - ۲02 

۲ 
اما به منظور تکمیل حل عددی نیاز به انتگرال- 
گیری زمانی معادلات (۱۲) و(۱۳) می‌باشد. که در 
تحقیق حاضر از روش نیومارک عمومی 022 برای 
گسسته زمانی پارامتر جابجایی و روش 011 برای 
پارامتر فشار استفاده می‌شود [28]. در این روش فرض 
می شود متغیرها در زمان *] معلوم‌اند و باید در زمان 


۸۱ + ]1 < ۳۰1 محاسبه شوند. 
متا صًّ > بسا 


اشمتاش۵ + اشنا + متا < وبملاً 

۱( ت"» مل و0 < + تارف < + شیرتا + وتا < ببسلاً 

+ ح < وبوق 

)۱۵ اشم6۵۲ + شم + «ط < دبوط 

بااعمال روش نموی(10076۳060121) تنها 

مجهولات این معادلات ما2 و ,۸ خواهند بود. پس 

از گسسته سازی زمانی» معادلات به صورت زیر در 
مایت 

2-0 ۳۳ - ۵ :0۳*16 - (۳*1) ۲ + متش۷ 

- رد + ۷۱۳۳۵۵۵ + رتاشاش,01(۳۲86) 

11 (۱10 

پارامترهای »0 و و8 در محدوده ۰ تا ۱ قرار دارند 

و برای پایداری بی‌قید و شرط روش حل باید 

۰ < 8 و 0.5 < 8 باشند [29], در نهایت شکل 

ماتریسی معادلات کوپله خطی به صورت زیر حاصل 


می شود. 
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۳۸ بکارگیری روش المان‌محدود در تحلیل همبسته سدهای تحاکی و برآورد فشا رآب‌متفذی 


مبملها اش ویرل] اوق + ربم0] 

مبو] ,۵8 +5 ببو 8,۵0 

1 1 
۳ - رم 

بطوریکه 

وت فمیاا] تیا وی( 1۳۳ > ۳۳ 

ی 
ص۱۱ ] یز ۲اعصری ۳۳۰۱ 

اش( ولا وبم[0] 
که ال طول گام زمانی است. دستگاه معادلات 
غیرخطی (۱۷) می تواند با استفاده از یک روش مناسب 
همچون روش تکرار نیوتون حل گردد. البته قابل ذکر 
است با فرض رفتار خطی برای خاک در یک تکرار 
جواب حاصل می‌شود. با حذف عبارات مربوط به 
شتاب (نیروهای اینرسی) از معادلات (۱۲) و (۱۳)» 
همچنین صرفنظر کردن از تراکم پذیری سیال و ذرات 
جامد. فرم شبه‌استاتیکی معادلات به دست می‌اید که 
اغلب در مسائل مربوط به تحکیم و محاسبه نشست 


کاربرد دارد. 
۸ () < (۱0[)5 - (0)[م] 
(۱٩)‏ (ر6) - (۴)[,] + ان 0[۳] 


در این حالت ,۲ وزن لا به خاکریز افزوده شده در 
هر مرحله و 0 < و است. بعد از گسسته‌سازی و 
نوشتن معادلات به فرم نموی. شکل ماتریسی معادلات 
شبه استاتیکی به صورت معادله (۲۰) خواهد بود [30]. 


نت تک ۱ ۳۳ كِِ ۱ تبو[0] مبول] 
(۸۵۱]۱)0 (۵)) [یبول] اش ببی"[0] 


ره 

لازم به ذکر است برای حفظ حالت تقارن ماتریس 
ضرایب. خط دوم معادلات در یک منفی ضرب شده 
است. [۵) بردار افزايش تغییرمکان گرهی و (۵0) 
اضافه فشار حفره‌ای ایجاد شده است. خط اول معادله 
(۲۰) تغییر در نیروهای داخلی (نیروهای ناشی از تنش 
موثر و فشار آب حفره ای) را به تغییر در نیروهای 
خارجی در گام افزایش زمان ۵۶ مربوط می‌کند. 
ماتریس سختی [ «ع] مرتبط کننده جابجایی‌های گرهی 
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به نیروهای گرهی است و با استفاده از رابطه تنش- 
کرنش اسکلت خاک حاصل می‌شود. بنابراین 
[۵][م] بیانگر نیروهای ناشی از تغییر تنش موثر 
می‌باشد. ترانهاده ماتریس اندرکنش» [0] مرتبط کننده 
فشارحفره‌ای در گرهها 
اختلاف »وبع(۳) و :(۲) برابر تغییر فشار آب حفره‌ای 
در مدت زمان ۸6 است. بنابراین اختلاف »ربم(۳) *[0] 


به نیروهای گرهی است و 


و ,(0[*)0] نشان دهنده نیروی ناشی از تغییر فشار آب 
حفره ای است. خحط دوم معادله (۲۰) بیانگر معادله 
تزازشی است. لازم هه دک ایک برای. کنست‌بازق 
معادله تراوش گذرا در حوزه زمان از تفاضل محدود 
استفاده شده است. اگر متغیر دورنیابی3 برابر یک 
باشد» روش پیشرو زمانی کاملا ضمنی بوده و بی قید و 
شرط پایدای است ۵1 القه عا باقن نظن گرفن 
5 < 9 معادلات را به صورت ضمنی حل خواهیم 
کرد. در مقاله حاضر برنامه‌ای براساس روش المان 
محدود در محیط فرترن توسط نگارنده توسعه داده شده 
است تا در طی مراحل ساخت بدنه سد. با حل معادلات 
حاکم و محاسبه مقدار افزايش تغییرمکان‌ها و اضافه 
فشار آب حفره‌ای» تغییرمکان‌ها و اضافه فشار آب 
حفره‌ای ایجاد شده در بدنه سدخاکی و پی آن را در هر 
گام زمانی محاسبه کند. بدین منظور ابتدا بوسیله 
سابروتینی داخل برنامه. بدنه و پی سد مانند شکل (۲) 
مش‌بندی شده و پس از اعمال شرایط مرزی و با 
استفاده از المان‌های چهار گرهی برای فشار و۸ گرهی 
برای تغییر مکان. ماتریس ضرایب و نیروها تشکیل شده 
و دستگاه معادلات ماتریسی حل می‌شود. خروجی 
برنامه شامل تغییرمکان‌ها و اضافه فشار منفذی ایجاد 


شده در هر گام زمانی اننتت: 


مشخصات سد مورد بررسی 
سد درونگر سد خاکی با هسته رسی قائم می‌باشد که 
در ۳۵ کیلومتری شهرستان درگز واقع در خراسان 
رضوی احداث شده است. ارتفاع سد از پی ۵۵ متر و از 
بستر ۶۰ متر می باشد[ ۳۲]. طول تاج سد 1۲۸ متر بوده 
و نسبت طول به ارتفاع سد ۱۰/۷ است. 


سال سی و یکم. شماره دو ۱۳۹۷ 


فرهود کلاته- فریده حسین نژاد 


تراز تاج سد ۹۰۲متر و تراز نرمال آب ۸۹٩‏ متر 
می‌باشد. خصوصیات و ویژگی‌های مصالح بدنه سد و 
پی آن در جدول (۱) ارائه شده است [۳۳و۳۶]. مقطع ۷۵ 
به عنوان بلندترین مقطع شامل ابزار دفیق برای تحلیل‌ها 
انتخاب شده است. شکل (۱) نشان دهنده ناحیه‌بندی 
پی و بدنه مقطع ۷۵ به همراه محل ابزاردقیق نصب شده 
دز آن می‌باشد. در این مقطع سه دستگاه پیزومتر 
الکتریکی در پی ((48 سه دستگاه پیزومتر الکتریکی در 
خاکریز (۳) و هفت دستگاه پیزومتر لوله ایستاده (5۳) 
نصب شده بود که متاسفانه پیزومترهای لوله ایستاده 
دچار نقص فنی شده و از کار افتاده بردند. 


۳۹ 


فرضیات محاسبه و پارامترهای تحلیل 

به علت طویل بودن تاج سد. محاسبات بر اساس حالت 
کرنش صفحه‌ای برای مرتفع‌ترین مقطع سد انجام شده 
و برای ساده سازی محاسبات. از تاثیرات سه بعدی 
بودن جسم سد صرفنظر شده است. مشاهدات علمی 
نشان داده‌اند که تحلیل دوبعدی تقریب بسیار خوبی از 
واقعیت می‌باشد. تنها در مورد سدهایی که در دره های 
تنگ با احتمال پدیده قوس‌شدگی ساخته می‌شوند. 
توصیه به استفاده از تحلیل سه بعدی اجزای محدود 
هی گر فة [35]. 


جدول ۱ مشخصات مصالح سد درونگر [22] 


نوع مصالح هسته رسی ‏ | پوسته شن و ماسه‌ای | پی آبرفتی 
شرایط زهکشی زهکشی نشده زهکشی شده زهکشی شده 
مدول الاستیسیته (1602) ۱۳۳۰۰ 3 ۱۷۳۰۰ 
ریت باس ۵ ۲/ ۳ 
نفوذپذیری قائم (0/9) ۳ ۳ ۳ 
نفوذپذیری افقی(0/9) ۳ ۳ ۳ 
مدول بالک ذرات جامد خاک(۵) ۲۰۱۵9 ۲.۱۵9 ۲۰۵9 
مدول بالک سیال(۵) ۱۰۵20 ۱۰۰20 ۱۰۰20 
چگالی ذرات جامد (12/۳3) ۷.۰ ۲۱ ۱۹۰۰ 
چگالی سیال (1>2/73) ۱۰ ۱ ۱ 
1 
۱ 


1 
‌ 


شکل ۱ جانمایی ابزار دقیق در مقطع ۸-۷ سد درونگر 


11 


سک 


260.93 


66۷۸۲۱۵۷ ۷ ۵۵ 


سک 95 


شکل ۲ مش بندی به کار رفته برای سد وپی زیر آن 


سال سی و یکم. شمارهٌ دو» ۱۳۹۷ 
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۳۰ بکارگیری روش المان‌محدود در تحلیل همبسته سدهای تحاکی و برآورد فشا رآب‌مفذی 


مطالعات انجام شده بر روی سدهای خاکی نشان داده 
است که با تو حه به نست 1/۳۲ (طول سد به ارتفا 

ی یه نمی زرم 
آن» باید درباره آبعدی یا ۳بعدی مدل نمودن سد 
از آن است که اگر نسبت طول به ارتفاع در سدهای 
اک مرن ایا ات فان ارم وی 
سد را به صورت ۲ بعدی مدل نمود [36]. این کار 
گردد. در سد مورد بررسی ما نیز همانطور که گفته شد 
مدل دوبعدی بهره برد. به طور کلی در تحلیل و طراحی 
سدهای خاکی. فرض مدلسازی کل سازه در یک 
قابلیت ساخت مرحله‌ای بهره جحست [۱۸]. تحقیفات 
نشان می‌دهد که مدل‌سازی در ۸ الی ۱۰ لایه می‌تواند 
به طور قابل قبولی جواب‌های منطقی و منطبق بر 
واقعیت ارائه دهد [37]. در این بررسی ساخت سد در 
۰ مرحله با بازه‌های زمانی مختلف و در مدت زمان 
کل ۸ ماه مدل شده است (حدول(۲)). برای این منظور 
مدل با شبکه ثابت در نظر گرفته شده به این صورت که 
وکا تاش مه وی از 
ساخت. وزن واقعی لابه اعمال شده است. آبگیری 
مخزن نیز ۱۰ ماه پس از پایان ساخت بدنه صورت 
گرفته شده و سطح آب طی دو مرحله تا تراز ۸۸۶۷۶ 
بالا برده ۷ بدنه اصلی سد در مدل. شامل هسته 
شده است. لازم به ذکر است در مورد مشخصات 
مصالحء پراساس محد‌وده نتایج آزمایشگاهی. نامناست- 
ترین ترکیب از دیدگاه ایجاد فشارهای منفذی در پی و 


بدنه در نظر گرفته شده انستا: 
شرایط مرزی 


به دلیل واقع شدن بخش تحتانی پی بر روی بستر 
سنگی» تغییرمکان های افقی و قائم مربوطه مقید شده- 
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اند. در دیواره‌های قائم پی تغییرمکان‌های افقی مقید 
شده‌اند ولی تغییرمکان‌های قائم آزادند. اضافه فشار آب 
حفره‌ای در دیواره‌های پوسته صفر منظور گردیده است. 
همچنین همزمان با آبگیری نیروی هیدرواستاتیکی فشار 
آب بر شیب بالادست و سطح افقی پی زیر مخزن وارد 
می‌گردد. نیروی دیگری که حین آبگیری بر نواحی زیر 
سطح آزاد آب وارد می‌گردد نیروی شناوری می‌باشد. 
سطح آب قیل. از آبگیزی, در تراز ۸1٩‏ فرار داشتة و 
پس از آبگیری تراز آب در بالادست به ۸۸۶۰۷۶ رسیده 
بود. 


جدول ۲ مراحل ساخت بدنه سد 
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تحلیل نتایج وبحث 
با توجه به هدف اصلی پژوهش که مطالعه توسعه و 
میرایی فشار آب منفذی در هسته می‌باشد و به دلیل 
وابستگی آن به تغییرمکان‌های ایجاد شده در این بخش 
نتایج مربوط به اضافه فشار آب حفره‌ای و تغییرمکان- 
های بدست آمده از مدل‌های شبه‌استاتیکی (روش 
معمول در حل مسائل تحکیم و محاسبه نشست) و 
دینامیکی. همراه با داده‌های ثبت شده بوسیله ابزار دقیق 
ارائه شده‌اند. در مقطع 2 از سد مورد بررسی. شش 
پیزومتر الکتریکی (۳8۳) در بدنه نصب شده است که 


شامل چهار پیزومتر در هسته و ۲ پیزومتر در پوسته می- 


سال سی و یکم شمارة دوه ۱۳۹۷ 


فرهود کلاته- فریده حسین ننژاد 


باشد. از پیزومترهای واقع در هسته نیز سه پیزومتر در 
راز رف فیگری دی ترا ۸۱۲ ققر تضیب شله‌اند: 
پیزومتر تراز ۸۲۲ (8۳7-7) در مرکز هسته واقع شده 
است و در تراز ۸۶۸ متر یکی از پیزومترها در مرکز 
هسته (8۳7-5) و دو مورد دیگر به طور متقارن در 
فاصله حدود ۱۲ متری از محور قائم مرکزی هسته قرار 
گرفته‌اند ( 12۳74 در بالادست هسته ۳7-6 در پایین 
دست هسته). که به دلیل طولائی‌تر بودن فاصله زهکشی 
در مرکز هسته نسبت به کناره‌ها انتظار می‌رفت مقدار 
اضافه فشار آب حفره‌ای در این بخش بیشتر از بقیه 
نقاط هم‌ترازش باشد ولی همانطور که در شکل (۳) 
مشاهده می‌شود بر خلاف آنچه که انتظار می‌رفت 
نمودار فشار آب منفذی مربوط به پیزومتر میانی پایین‌تر 
از دو پیزومتر کناری می‌باشد و ممکن است قرائت 
اولیه این پیزومتر اشتباه بوده باشد چون در ادامه مشاهده 
می‌شود روند موجود مطابق با روند نتایج حاصل از مدل 
المان محدود می‌باشد. 
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شکل ۳ نتایج حاصل از ابزار دقیق واقم در تراز ۸۶۸ متر. داعل 
هسته 
نتایج حاصل از مدل المان محدود می‌باشد. نتایج حاصل 
فلز ‌های, آنتاتکی. ور دیناسکن با قاایر لت 
نفوذیذیری در شکل ۶ نشان داده شده‌اند. همانطور که 
مشاهده می‌شود در هردو تراز مورد بررسی با کاهش 
نفوذپذیری مصالحء مقدار فشار منفذی ایجاد شده در 
هسته افزایش یافته و به مقادیر مشاهداتی ابزار دقیق 


۰ ۳ 1 هر ۱ 
نزدیک‌تر می‌شود به طوری که در نفوذپذیری ‏ ۱۰ 


سال سی و یکم. شمارهٌ دو» ۱۳۹۷ 


۳۱ 


نتایج مطابقت خوبی با داده‌های ابزار دقیق دارند. 
همانطور که ملاحظه می‌شود در تمامی حالات نتایج 
حاصل از تحلیل دینامیکی بیشتر از تحلیل استاتیکی 
بوده و اختلاف موجود بین نتایج این دو نوع تحلیل با 
کاهش نفوذپذیری افزایش می‌یابد. علت را اینطور می- 
تون نیاق کرد کهدن سالت شبه‌استایکن از تار شفاب 
ناشی از سرعت بارگذاری صرفنظر می‌شود در حالی که 
در ساخت لایه‌ای بدنه فرض بر این است که کل وزن 
لایه در ابتدای بازه زمانی مربوط به ساخت آن لایه 
اعمال شده و در باقی بازه زمانی زهکشی صورت می- 
گیرد و همین امر موجب ایجاد شتاب شده و نتایج را 
تحت تاثیر قرار خواهد داد. مطابق انچه که در شکل 
(۶) مشاهده می‌شود برای هر کدام از ترازها اعتلاف 
بین نتایج دو نوع تحلیل در زمان‌های اولیه ناچیز بوده و 
با گذشت زمان یعنی افزايش ضخامت خاکریزی بیشتر 
می‌شود و پس از اتمام ساخت بدنه به دلیل زهکشی 
اندکی کاهش یافته و در ادامه ثابت می‌ماند. علت آن را 
این طور می‌توان توجیه کرد که وقتی ضخامت 
خاکریزی کم است مسافت زهکشی قائم در تراز مورد 
بررسی کم بوده و علاوه بر زهکشی افقی» زهکشی قائم 
نیز اتفاق می‌افتد ولی با افزایش ضخامت لایه‌ها امکان 
زهکشی قائم برای لایه‌های زیرین کاهش يافته و اضافه 
فشار منفذی ایجاد شده تقریبا فقط از طریق زهکشی 
افقی صورت می‌گیرد. در شکل‌های (۵) و (1) توزیع 
کنتوری نتایج حاصل از هر دو نوع تحلیل برای 
تغییرمکان‌ها و فشار حفره‌ای ایجاد شده پس از اتمام 
ساخت بدنه و بالافاصله پس از آبگیری نشان داده شده 
است همانطور که ملاحظه می‌شود الگوی توزیع 
تغییرمکان‌ها و فشار حفره‌ای در هر دو نوع تحلیل 
یکسان بوده وتنها تفاوتشان بیشتر بودن مقادیر حاصل از 
یا شتاستکی سب آسفاشکی اس من فتاه 
قرف ریق دای فشار یی تاه شقن بات 
در محل اتصال هسته به پی می‌باشد که دارای بیشترین 
مسافت زهکشی است. در یک سوم بالایی هسته و در 
نزدیکی فیلترها به دلیل امکان زهکشی بیشتر. اضافه 


فشار منفذی ایجاد شده ناچیز است. همچنین به دلیل 


نشریهٌ مهندسی عمران فردوسی 


۳۲ بکارگیری روش المان‌محدود در تحلیل همبسته سدهای نعاکی و برآورد فشا رآب‌مفذی 


نفوذپذیری بالای مصالح پوسته هیچگونه اضافه فشاری 
در آن مشاهده نمی‌شود. مقدار تغییرمکان‌های افقی که 
متاثر از شیب دیواره‌ها می‌باشند به دلیل شیب نسبتا 
مشابه شیروانی بالادست و پایین‌دست کاملا متقارن بوده 
و همزمان با آبگیری اندکی به سمت پایین دست متمایل 
می‌شوند. تغییرمکان های قائم نیز پس از آبگیری به 
دلیل نیروی شناوری وارد بر مصالح تقلیل می‌یابند. پس 
از آبگیری مقادیر اضافه فشار حفره‌ای نیز کاهش می- 
یابد. به منظور بررسی زمانمند تغییرمکان‌ها و فشار آب 
حفره‌ای حاصل از هر دو نوع تحلیل و تاثیر نسبت 
نفوذپذیری افقی به قائم در نتایج تاریخچه زمانی 
تغییرمکان‌های قائم و فشارحفره‌ای در دو تراز ۸۶۸ و 
۲ متر برای نسبت‌های مختلف نفوذپذیری افقی به 
قائم به همراه نتایج ابزار دقیق به ترتیب در شکل‌های 
( و( نشان. داخه اند متانظور. که شام 
می‌شود در مرکز هسته با افزايش تراز خاکریز تاثیر 
نسبت نفوذپذیری افقی به قائم مشهودتر است. به 
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طوریکه حساسیت نتایج به نسبت نفوذپذیری افقی به 
قائم در تراز ۸۸۲ بیشتر از تراز ۸۶۸ بوده است. علت را 
می‌توان کم شدن ضخامت هسته با افزايش تراز و در 
نتیجه کم‌تر شدن مسافت زهکشی افقی در ترازهای 
بالاتر دانست. همچنین در یک تراز مشخص نیز به دلیل 
کاهش مسافت زهکشی, این تاثیر با پیشروی به سمت 
دیواره‌های هسته بیشتر می‌شود. به طوریکه در تراز ۸۶۸ 
تغییر نتایج در بالادست و پایین دست بیشتر از مرکز 
بوده است. لازم به ذکر است نتایجی که در آزمایشگاه 
برای نفوذپذیری خاک بدست می‌آید. به نفوذپذیری 
قائم خاک نسبت داده می شود و نفوذپذیری افقی بسته 
به شرایط موجود بین ۲ تا ۱۵ برابر این مقدار در نظر 
گرفته می‌شود [۱۸]. 
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شکل ۶ مقایسه اضافه فشار آب حفره ای بدست آمده در دو حالت تحلیل شبه استاتیکی و دینامیکی با داده های ابزار دقیق 


و تراز ۸۶۸ متر(بالادست هسته), 0. تراز ۸۶۸ متر(مرکز هسته). ». تراز ۸۶۸ متر(پایین‌دست هسته) 0. تراز ۸۱۲ متر(مرکز هسته) 


نشرية مهندسی عمران فردوسی 


سال سی و یکم. شمارهُ دو. ۱۳۹۷ 


فرهود کلاته- فریده حسین نژاد ۳ 


شکل ۵ نمودار توزیع کنتوری تغییرمکان‌ها و اضافه فشار حفره‌ای ایجاد شده در مقطع عرضی سد در انتهای ساخت برای دو حالت تحلیل 
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شکل 1 نمودار توزیع کنتبوری تغییرمکان‌ها و اضافه فشار حفره‌ای ایجاد شده در مقطع عرضی سد بعد از آبگیری برای دو حالت تحلیل شبه 
استاتیکی(سمت چپ) و دینامیکی(سمت راست) 


سال سی و یکم. شماره دو. ۱۳۹۷ نشریهٌ مهندسی عمران فردوسی 
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بکارگیری روش المان‌محدود در تحلیل همبسته سدهای خحاکی و برآورد فشا رآب‌مفذی 
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شکل ۷ نمودار نشست‌های بدست آمده از مدل المان محدود و داده های ابزار دقیق در ترازهای ۸۶۸ (012161) و ۸۲ (5عاقام) متر برای دو 


جدول ۳ درصد اختلاف نتایج تحلیل و داده‌های ابزار دقیق در دو مرحله پایان ساختمان و آبگیری 
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در اینجا نیز نفوذپذیری بدست آمده در آزمایشگاه 
به عنوان ورودی اولیه برای شروع تحلیل‌ها به کار رفت 
و در ادامه به منظور مطابقت بهتر نتایج حاصل از مدل 
با نتایج ابزار دقیق با توجه به نتایج حاصل از نسبت‌های 
مختلف نفوذپذیری قائم به افقی, نفوذپذیری افقی 
حدود ۲ برابر نفوذپذیری قائم خاک در نظر گرفته شد. 
لذا مقدار نهایی بدست آمده برای نفوذپذیری افقی و 
قائم مصالح هسته در تحلیل‌های برگشتی به ترتیب 
عبارت اند از: 
6/960۰ 107۲ <1 رک ,(۳/566) 10711 <2 ع بک1 


لبته برای تراز ۸6۸ در مرکز هسته تغییر 


نشریهةُ مهندسی عمران فردوسی 


تراز ارتفاعی 
ِ_ِ ی 
10۵ نوع تحلیل 
53 لها 
۰ 2292010 
81 متصحصرل 
مها 
2-2211" 
متصعصرل 


نفوذپذیری افقی تاثیر چندانی بر نمودار نداشت ولی در 
تراز ۸۳۲ انتخاب ,26 نتایج نسبتاً مطلوب‌تری 
(شکل ( و (۸). با 
استفاده از ضرایب نفوذپذیری فوق مقادیر فشارهای آب 
حفره‌ای اندازه‌گیری شده و محاسبه شده برای مدل 
المان محدود هماهنگی بسیار خوبی با هم داشتند. به 
طوریکه در زمان پایان ساخت مقادیر حداکثر فشار آب 
مفره‌ای سفاصل از قحلیل شبه استخاتیکن حدود ۱۳ ۶ ۱ 
درصد و نتایج حاصل از تحلیل دینامیکی حدود ۳۰ تا 


۵ درصد بیشتر از داده‌های ابزار دقیق بودند (جدول۲). 


تسه .تیه عتالت‌ها :بویت 


سال سی و یکم. شمارهُ دو. ۱۳۹۷ 


فرهود کلاته- فریده حسین ناد 
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0 ( چ ۴ 0 + + 50000 
۳ +7 40000 وت 3 غِ 
ص۳5 ط 70 3 هه 3 
مومع + 30000 1 توت ج مومود ع 
2۷ تب ورو [ 1 ود ۵ 260 1 20000 
تست ۳۶ [ ی سس 850 7 + 100,00 
65 ۲ ومع 1 ات هسب 00 
+ 00 2 9 1 ى‌ِ 
۳ ست. 860 7۴ ‌ ۳ فن 2 2 ۳ هِ 2 دِ 
این اه 2 1 1 اععا ۱ 8 5 5 ل 5 43 
وج 8 مه و و 8 چ و ٩‏ اه دص 5 
سست گ 5 ج چ چ شست چ جع ی ۲ 3 
تفت چ ۶ 5 5 2۷ 9 2 
۳ ی 0 7 910 > 80000 
910 7 0000 7 
موه 67-5 سس 900 تن 
سس و نت ۱ 60000 
500,00 تست 
ادوس ووع ٍ 907۳ + 50000 و 
+ 40000 
خ اسب + 40000 
اس ووو ا اس 0و8 : ِ 
1 ۱ لا اس 70و ۶ ومممو 1 
سس ۱ 50| سس 3 
0 ۱ 7 200.00 + 20000 عه 
: 3 
سس م6 ۱ 10000 سس 50و 10000 
سس 850 ‌ ۱ سس 850 ۳ ست ححد تج سچ. 000 
لت ۳ لت ۱ 1 هِ 2 1 ۳۹ 
یت سر ۸ ۰ ۰ ۰ 3 ۳ 
2۷ 5 3 4 3 3 
910 0000 
1 20000 
است 7-6] سس 900 ۱ 
76 ووو 1 ت__ 
اس وچ جاعتاس 990 نت 
2200 سس 80 تشخ : 


ا5-اس 70و + 


۳۱۳۰ )۱( 


سس 860 


سس 85 3 
سس 850 


(هها) ۳ 
1 
3 


4 ۰ 000 
شنت ۲ 9 
وج 840 ۸8 سس و۵4 27 3 3 10000 
8 اسب ۵ 4 "_" 4 ۵ 
۷ ِ ك ۳ 
05 50000 
910 1 
1 دوعس 9990 [ 50000 
7وست. 1 
1 ووو + 
ین ۱ (0-بواسب- 1 40000 
نات 1 0 + 1( 
0 1 نایبت 1 ت 
سای .۱۰ چ 5و + ۲ 
1 غِ << ووع 1 20000 
880-5۷ 1 ِ 1 و 
1 200 وزو + ا تفت 
اس ووو ‏ ِ 
ون 1 ۳ ۰ + 000 
1 3 بچبت 
عت ۰ اعبع۱؛ 3 
ا۵ع . 850 7 1 65 چٌ 4 ۸ ۸ ِ 1 
۳ 4 یبا 3 1 ۱21 1 1 10000 
سب ۳ س__ مه 5 3 ِ ِ ِ 1 
250 1 ۳ 5 5 2 3 
3 ِ 5 . 5 ۶ 3 


شکل ۸ نمودار فشار آب حفره ای بدست آمده از مدل المان محدود و داده های ابزار دقیق در ترازهای ۸۶۸ و ۸۱۲ 


نتیجه گیری و جمع‌بندی باشد و نتایج حاصل حاکی از این است که: 
هدف از این مطالعه. توسعه و حل معادلات 80۶ با ۱. کلیه نمودارهای فشار آب منفذی و تغییرمکان روند 
استفاده از روش المان محدود جهت محاسبه فشار افزایشی داشتند که با توجه به روند خاکریزی و 
منفذی وتغییرمکان‌ها در هسته رسی سدهای خاکی می-. افزایش سربار منطقی می باشد. 


نشريهة مهندسی عمران فردوسی 


سال سی و یکم. شمارهُ دو. ۱۳۹۷ 


۳۹ بکارگیری روش المان‌محدود در تحلیل همبسته سدهای تحاکی و برآورد فشا رآب‌مفذی 


۲ مقدار نفوذ پذیری افقی در تراز های بالاتر نتایج را 
بیشتر تحت تاثیر قرار می‌دهد. مقدار نهایی که در 
تحلیل‌های برگشتی برای نفوذپذیری افقی و قائم 
مصالح هسته بدست امده به ترتیب عبارتند از: 
(۳/560) 1011 1 > ,(560/) 101۲ <2 ,1 
۳ میزان فشار آب حفره‌ای ایجاد شده و همچنین 
سرعت افزایش آنها در هنگام خاکریزی در نقاط 
شاف اه ویس به: ظراملی تفانتد قراز نشب 
پیزومترهاء فاصله پیزومترها تا فیلترها و... بستگی دارد 
به طوریکه در یک تراز مشخص. بیشترین مقدار فشار 
آب منفذی در بخش میانی هسته بوده و به سمت 
کناره‌هاء از این مقدار کاسته می‌شد که کاملا منطبق بر 
واقعیت می‌باشد. 
۶ نتایج حاصل از تحلیل دینامیکی به طور متوسط ۲۰ 
درصد بیش از مقادیر بدست امده در تحلیل شبه 
استاتیکی بودند. 
۵. مقادیر حداکثر فشارهای آب حفره ای ایجاد شده در 
زهانتبابای اه هیک اش اتکی 9۱۳ ۱ 
درصد با داده ابزار دقیق اختلاف داشتند و در تحلیل 
دینامیکی نتایج حاصل از مدل ۲۰ تا ۳۵ درصد بیشتر از 
داده های ابزار دقیق بودند. 
1 در صورت بالا بودن سرعت اجرای خاکریز استفاده 
از معادلات دینامیکی ضروری است و تنها در صورت 
سرعت پایین اجرای خاکریز نتایج حاصل از تحلیل 
ساخت مرحله‌ای سد با معادلات شبه استاتیکی قابل 
قبول خواهند بود. لذا می‌توان با در نظر گرفتن ساخحت 
مرحله‌ای سد در تعداد لایه‌های بیشتر و لحاظ کردن 
زمان ساخحت نسبتا طولانی مراحل ساخت سد را با 
مراچع 


معادلات شبه‌استاتیکی تحلیل نمود. 

یک مدول بالک دانه های جامد 

لا شتاب اسکلت جامد 

جابجایی اسکلت جامد 
سرعت نسبی جریان سیال نسبت به 
فاز حامد 

بو شتاب جریان سیال نسبت به فاز جامد 
زط . نیروی حجمی وارد بر واحد جرم 
وتات 

وی نفوذپذیری دینامیکی 

2 بل 

و فشار آب حفره ای 

1 چکالی حفرات 

2 چگالی دانه های خاک 

5 

0 چکالی اب 

۷ وزن مخصوص آب 

ی تنش کل 

ع کرنش 

مدول بالک سیال 


0 ۲۱۷۱۵0۵۵۵۵ :0۲۵۵۵۵(۰)) مرچ‌نز طاتوظ ومصیما۱ مطا ۵۶ وتووامصق کصمجصما)اعو" رب رتاملت2112 .۷ رفقئلات 


و0۰11 و30 ۷۵۱ رکه امه وصامع اوه دیول و فصتماندم۱۷ متاع00ع0 مه عصتام۱۵0 اهمزمصصیا 
2008(۰) .3074-3081 .00 


۲ شرکت مهندسان مشاور آشناب "راهنمای تحلیل و طراحی لرزه‌ای سدهای خاکی و سنگریز؛معاونت نظارت راهبردی, امور نظام فنی. 


نشریه شمارهء ۰1۲ (۱۳۹۲). 


6 0۱۲۱۱۵ ۳ظ1۱(۵1018۵)10 صرح وعماه امممص‌اموماط ۵۶ وتورامصم ۲۷۷۵۵02 .ل بر و۷۷ .0 مرطعتامن 
٩۲0, 4, 00. 529-549, )1967(,‏ و93 ,۷۵1 وطمکذ۲ ععه0م۲090 090 1۷/6۵۲۵۳165 5011 


نشرية مهندسی عمران فردوسی 


سال سی و یکم. شماره دو ۱۳۹۷ 


فرهود کلاته- فریده حسین نژاد 


۳۷ 


۶ زمردیان س.م.ع. و چوچی. ح. «تحلیل عددی رفتار سدهای خاکی- سنگریزه‌ای حین ساخت و اولین آبگیری(مطالیة موردی سد 


مسجدسلیمان)». مجله علوم و فنون کشاورزی و منابع طبیعی, علوم آب و خحاک:7۲(۱).صص. ۰۲۲۹-۲۶۲ (۱۳۹۱). 
کصمحصلصح ماحصظ هن و ۵۶ ممتامنعومن مط) مصتلام۰۱۷]۵0 رت و۲۱ ععم؟ رب رعطه آ و.۲ و20 ۲۲0۵ ,.۵ مصفط) 
(2014) ,93-102 00۰ و0.1ظ و18 ۷۵۱ و19 ۱66۲و از 0 [۵یاهل نی ر"صصور1 
هلدابم ماجنا امممصتممدظ مه فممتامن۳۳۵0 فتوراممم اهمزمصصتان مصتنهم‌صمن رب رلهاکض۸ 
2009(۰) 50-62 .00 با ,0 و2 ۷۵ ,56۱6۳065 ۱۵0ا۵۵ ]0 0۲۵1 مزع صصجنا نع کصمحصصه‌ماهر۴ 
جح گم بناامدنه مجح تصمصه‌ااامهو ۵۶ وتوراممه مه مصتممانمم۱ بلط مرطعقصهآ. رفظ رلتامجطم]۲]-آظ 
6 00۱4 ب«و010ع 01۳۳9۱۳۱۵۵۲۱۸9 8۰ع۱ابظ و0صنامن 5011 مه ومعهاه صا عقوم ون اصممص‌اصو‌طماصر۴ 
(2001) ,259-2864 .00 و04 و59 ۷۵1 ۲۱۱۷۱۲۵۳۱۵6 
۷۵۱16 #مع۲ جم حصونا عاممعطاتجظ ۵۶ ممتاقططمگونا ووعتاو ۵۴ فلورفمد تحمجصماط مالصا۳۳۲ ر.۲ بت 
(2014) ,1852-1855 .00 ,580-583 وکاهز ۱۷/۵6۲ ۵8 کعن۷۲6۵۳۵۳] 0عز ام 
وبا ۸ وکاممدوملهو بلط رئتهلمصاهتهطه ری ۱۷ وتقاموظ ویک با مصهتمطمهها ت۳۵ وی رن‌صفطف‌امرهرز 
۰ 296)) ومتامنتتاوومن) ۵ ۴80 مطا 2 ون عم ه ۵۶ مصازماتمم]۱ فصه فلو ولمم اهمزتممی۳۸ 
صمتامعوصع۳]. ولا و12 ۷۵ ,۱6۵9و ابش ۵ هیال م6۵۵0 و" (صعنا توروطج ۲۷ 
2014(۰) ,ع19661158عصظ! احعنصطه‌ما0ع0) 
۵ ها و -عووای وم مرها م9 کصمصص‌آصو‌مصطر ۵۶ امع۳ .۲ رلففصصاهه .۲ ات12 
2014(۰) ر369-375 .00 4 ٩0.‏ و5 ۷۵ 6۱۱۲۵ ]0۵ ]۳۵اه 
تمهت ممحتمو۲ ۵۶ عتوراممد ماو ری رتمعافلک ری رمصهطعتعظ ری رتصصهعمطت ۷ ,۲62266 
,۳6560۲6 6 ۱۱۵۷۵۵0۲ 1۳9۳6۵9 ]0 م0۵ ۱۱۵6۲۳۸۵۵۵۴۵ ,۱۷۵۱۵۵0۲ اصمجطماظ فانصا مصتونه حور 
(2014) ,654-658 .00 ,5668 - 2277 :199 و5 56و و3 ,۲۷۵1 
صتحتاک-عوم‌گ م1 وم ممتاقانه ۷ تطعتع۳۱ 0۵۴ ۳۶۶۵۵۲ ریک .۲ مط0مطقلظ یگ رتصطه‌عقدات رم مطعتصه ۱۷22 
0ص ۵ ۳۵۱00 ]0 آمتیا0ل 1۱0۱۵7 "عون کجممصصجطاص فناممممع۳۱۵)۵۲۵ 1۳ ۵۶ ولو رراهصض۸ 
۵ 0۱9۲۱۵۲( اوصمتاعصهاص مرن رووععع۸ مهم مه (فمتلمو0) 6345 -2231 :199 ,506765 ع]ز,1 
(2015) ,35-42 .00 و(94) 5 ,۵1 ۷ 6.187 /۵۱۵۵/۱۰۵۳۲9/5۲۰60/۱5/2015/04 01۱۷۷۱۸۰۵ 
هط بامآمطم‌صعنا مطا ۵۶ ۳۲۶۲۵۵۲ مضه مه ممنهمص۲۱(۵۲۵ روفعتاو ۶ ولوراممم م1 ,۷۷ رصمعط وبا رعط26 
2015(۰) ,5061617 [2عنصطمع)۵ع6) صهت‌صا رصصونا ات۲۵۵ 


10 


11۰ 


12 


13. 


3 زمردیان. س.ع.ع. و صاحب‌زاده. ک. (بررسی تعریف تعداد مراحل ساخت سدهای خاکی و سنگریزه‌ای با انجام آنالیز تحکیمی در 


مدل‌های کامپیوتری». مجموعه مقالات اولین کنگرهُ ملی مهندسی عمران دانشگاه صنعتی شریف. تهران» (۱۳۸۳). 


۵ پاکزاد.م. و حسینی. س. م. م. م.» «ارزیابی فشارهای آب منفذی در پی و بدنهٌ سد کرخه پس از آبگیری». مجموعه مقالات چهارمین 


کنفرانس سدسازی, تهران. (۱۳۷۹). 


سحسینی» س. م. م۰ م۰۰ موسوی خوانساری» م. و احمدی فرد. رضاء «ارزیابی فشار آب منفذی در هستةٌ سد کرخه هنگام ساختن و 


آبگیری هم‌زمان». مجلة مهندسی آب. سال دوم. بهار ۱۳۸۸ صص. ۳۷-۶۷ (۱۳۸۸). 


۷ نیرومند ح.» میرقاسمی. عَ. و پاکزاده م «رفتارنگاری سد کرخه در دوران ساخت بااستفاده از نتایج ابزاردقیق»» مجموعه مقالات 


چهارمین کنفرانس سدسازی, تهران (۱۳۷۹). 


۸ حبشی,ع.. بارانی.غ. و خبیری»و» «ارزیابی فشارهای آب منفذی در هستهٌ سد خاکی دوستی بااستفاده از آمار اندازه گیری‌شده و نتایج 


تحلیلی» مجموعه مقالات پنجمین کنگرء‌ملی مهندسی عمران. دانشگاه فردوسی مشهد. مشهد ایران» اردیبهشت. ۰۱۶-۱ (۱۳۸۹). 
مج ممتامنتتاعممن مصتتانا رللها۱ ۲ ما۲ [فها ۵۶ ممتماه ۲ اهگهو رتهگ1 7 جخ هر لمع ۷۲۰ رتعلع1۷]21 
٩۲0.101-85, )2005(,‏ رصق رصفعطه؟ ,]100 ۶و عهتاهع۱6 امنصهم 73۴ ,"هصنه‌ماهوا معما۹ اعیز۴ 


19۰ 


سال سی و یکم. شماره دو ۱۳۹۷ نشرية مهندسی عمران فردوسی 


۳۸ بکارگیری روش المان‌محدود در تحلیل همبسته سدهای تحاکی و برآورد فشا رآب‌مفذی 


0 سروش: ع. اخترپو ع.: «مطالعه و تحلیل عددی فشار آب حفره‌ایایجادشده در هسته سد سنگریزه‌ای مسجدسلیمان». مجموعه مقالات 


هفتمین کنگر؛‌ملی مهندسی عمران دانشگاه تربیت مدرس, تهران. (۱۳۸۵). 
مج جمتامصرم]عر۱ متصصتهصون فصع متاهاه مط صا 0مط)۱۷]۵ اصمجدرما مانصا۲ فظ] رفظ رهام .۴۱۷۷ و6715[ 
,(1998) ,۷۷۱۱6۷ :تعاععطمنط) ,60 2۳ ,2ن۱۸60 فبامتوظ و جمناعل‌تامووم) 
و12 ۷۵۱۰ رکعکت۳ظ ۱۵0ا م۸ ]0 امیاو0ل طمتاهلتامعجمی اقطمتعصمصصتل-ععتطا ۵۶ معط لهتمجهن" ره راماظ 
1941(۰) ,155-164 00۰ ,2 0۰ 
ام 066۵۲2112760 مطا به0۲ع۱ فیامعوظ 8ممتنهو گم تمزمطمظ متصهوونا رآ رتطملتا و0۵ و2م71001616771 
امعتاام( مه مع۱ ۳۵۱۸۱۵ ۵۲ آمیاول آمصمم 6 رصمتناامو. لجمزنممصیاط و105 مج صمتهآناطرن۲ 
.(1984) ,71-96 .00 وا ٩0.‏ و8 ,۷۵۱ رکع 660۳16۵۵0۳ 17 ۷۲6۱۳۱005[ 
۶ 2710۲ظ۵ متصتحصول 4صه متتهاو .۲ رتحطمنط بکلانا رلناعظ ,۷۲ رتماوهظ .۸۳۱ رصق و0 و7100161677162 
6 ۱۲۵6۵۵0۱9۶ ر"وصهاه9ازن عماج النا۳ .۲ .ممتانتامی مامتان ما طمهمتووع تقصمتاع1 2 وفلامو 
ب311-321 .00 ,(429)1877 ,56۱6۴1665 معا هه ۷۲۵۵۱۵۳۵۵6۵( رش 56۲۱۵5 .۵800 ]۵0 بچع506 ۲0۵ 
1990(۰) 
کمص م0۵ آمممتاهانام‌طمن ‏ رآ رتعطمتطه بفظ رتقااهتطهه ۲ رتمافوط ,۸۳۲۱ رصق و0 710010167712 
1999(۰) ,۷۷۱۱۵۷ :۷۵۲۱ ۱۱۵۷ ,"6صاعص ملمتمطا هط ما عمصمعام؟ آهزمع۵٩‏ ۳۷1 
متصصحصون ما عاع۱۷۲۵0 ی 0ععهظ اتمتاوهاظ ۵۶ مملاماصمصصمام‌صا ری .و تتمقا ریک یه تصصو‌ج2ظ ریک .۸ وله 
مق مضه .21۵0001۵2 مافج ۵۶ صموتنهم‌صصمن) . ه . :عون تمصع طاتعظ ‏ ۵ . ولورراحمض 
,(2004) ,385-410 .0 و3 ۷۵ بکه0)60ع6 ۵0 ۵۳۵6۲5 ]0 ۱۷۲۵06150۱۱۵ 


۱۷۲۵0۱2۵ عبامعوظ واماممرمگع۱ ۵۶ مصتاه۱۷]۵0 )ممحصماظ فص 460صم ۳۲ ریظ راقطعنطعه۲ یه رتمک 


2011(۰) ,5426-5441 .00 م11 .۱0 و35 ۷۵۱ ,1006۱۵ ۷۲۵۵۱6۳۴۳۵۵6۵1 ۸۳۵۵]۱60 ,ومع عاصا بتحتااهاتض 
0 ۲0۲ کمتا۵0عم۲ ممتاناآمو 0۲00معهاو ماماهای «ااممم0صممصنا و و0۵ ر2م1 2۱0۵1610 و۷ رعطع۲۲ 
,۱9۱۵6۲۱۵ ( ۷۵006 ۵۱۸۱۵۱۵۵ ۵۳ امدتیول ۲۳۸۱۵6۳۳۸۵۵۵08۵1 و عصجعامامنم متصصعصیول 1010 5011-20۶28 

(1998) ,1029-1052 .00 ,6 .0 ,43 ,۲۷۵1 
,100008 مطصهعهم ام ]رهاظ ر0مطاه]۱۷ عممعجما۴ فانصا م۳1 ول رتماره1 و0 71001016712 
2000(۰) ,561 .00 و1 :۲۷۵1 
همعط ,7110 , طمنانقه 4۳, ۱6300۲ کصمحرما مانصا مط مصتصهه۵ ۳۳ .1.۷ فطاتگنن ۲۸ طانصصه 
2004(۰) ,96-100 ۳0۰ 
متصطمم) و۲۵ هنصلع۷1 روتوعط]. ,قاط رعممتاه/۲۱۵۵۲ و۱۳۵ ۵۶ فلورافصض تممطورما۳ ماتطا۳ ر.ل.1۳ وتع04عظ 
(1998) ,83-87 .00 مهلطلع1 ۷ رعتادایمقاظ روتوم نصا متماه 0صه ماباتاعم] 
۳۲ مهنداسین,مشاور کاوش بی.مشهد: گزارفن ,رفتان نگازی مد (درونگر)::(۱۳۹۱): 
۳ مهندسین مشاور طوس آب. گزارش منابع قرضه و مصالح ساختمانی. مرحلدوم مطالعات سد مخزنی درونگر, (۱۳۸۳). 
۶ مهندسین مشاور طوس آب. گزارش فنی پی آبرفتی. بدنة سد و ابزاردقیق» مرحلدوم مطالعات سد مخزنی درونگر, (۱۳۸۳). 
۲6۵۷101۲ صوند وم ممتلهتاع‌لممن بومالج ۲۷ مطا ۵۶ وامع۳۲ رب رقصصههتنا له ,۲ رتتتتلنا۳ رب رهونک 
1993(۰) و95-106 00۰ مطمتاهنه۲۷ اجه ون وم م۵۲ ۲۷۷ امصمتعصتمام] 
و216 رک1۵0 عوم۱ ۵ ممکک۵۳ ۲۳۱۵۲۳۱۵۵۵08۵ و قصصونا غجمحصی‌صو‌مامر۴ ۵۶ وتورولفمد متاهاو" رثلنامن] 
.(1986) ,53 تاع۱11ظ 
.(1982) و2 ,۷۵۱ ۱۷۲۵۹۵0۷7 ر۲مکزامانام ۷ و منود متانهه ۳ رب رتیت 
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نشریهةٌ مهندسی عمران فردوسی سال سی و یکم شمارهُ دو ۱۳۹۷ 


